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Abstract 

We report a collection of lowest-energy structures of hydrocarbon molecules C m H n (m=l-18; 
n=0-2m+2). The structures are examined within a wide hydrogen chemical potential range. The 
genetic algorithm combined with Brenner's empirical potential is applied for the search. The 
resultant low-energy structures are further studied by ab initio quantum chemical calculations. The 
lowest-energy structures are presented with several additional low-energy structures for comparison. 
The results are expected to provide useful information for some unresolved astronomical spectra 
and the nucleation of growth of nano-diamond him. 
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INTRODUCTION 



Hydrocarbon molecules, especially polycyclic aromatic hydrocarbons (PAHs), have re- 
ceived much attention for the broad astronomical interests. PAHs were proposed to be the 
possible carriers of the unidentified infrared (IR) emission bands[l|, the interstellar ultra- 
violet (UV) extinction curve 

as 

and the diffuse interstellar bands 0-0]. Though quite a 
few achievements have been made in recent years [3, Q], the PAH hypotheses are not fully 
addressed. Possible different charged states of PAHs may increase its complexity, however, 
one intrinsic difficulty could be stemmed from the very complicated phase space of hydrocar- 
bons due to the enormous bonding ability of carbon. Furthermore, the phase space is also 
dependent of the hydrogen chemical potential, as the observed IR spectra may be different 
in different interstellar environments [8|. 

In this paper, we focus on providing a collection of lowest-energy structures for C m H n 
(m=l-18; n=0-2m+2) in a wide hydrogen chemical potential range with the aid of a series 
of unbiased global searches. Our results are expected to be helpful for addressing the above 
PAH hypotheses. Moreover, since the most abundant reactive elements in the cosmos are 
carbon and hydrogen, our database should be helpful in elucidating the evolution of carbon 
from its birthsite in circumstellar shells through interstellar medium, which may be related to 



astrobiology i lit in lately [8]. Useful information for a better understanding of the mechanism 
of^nucleation and growth for nano-diamond films may also be extracted from our results^, 



10|. 



COMPUTATIONAL METHODS 

T he g lobal structure optimization is based on our early developed genetic algorithm 
(GA) [HI]. An illustrative online version of the code for the optimization of a two-dimensional 



map has been published in nanohub 12]. Each GA run is started with a randomly gener- 



ated pool with typica 
empirical potential 



ly 100 structures which are relaxed to local minimum with Brenner's 
131 ] . The size of pool may vary with different investigated systems. The 
random generation of the pool may be guided by some physical intuition. The evolution of 
the pool is realized by performing mating operations. Each operation includes randomly se- 
lecting two structures (parents) from the pool, cutting them with a common plane randomly 
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selected, then combining the opposite parts relative to the cutting plane from the parents 
to create a new structure (child). The child structure is relaxed within Brenner's modeljisl]. 
The decision of replacing the highest-energy structure in the pool with the child structure 
is based on energy criterion / defined as 

f = (E - n H /j, H )/n c (1) 

where E is the total energy of the cluster with nn hydrogen atoms and tlq carbon atoms, nc 
is fixed for each GA run, while n# may be varied in a physically reasonable range. /z# is the 
hydrogen chemical potential, which may be fixed to some particular value(s) [14] or a range 
of interest. In this work, we are interested in a wide hydrogen chemical potential range, 
which is taken as (-6eV, 0) in Brenner's model, i.e., the lowest (highest )-hydrogen chemical 
potential value guarantees that the optimum molecules are pure carbon clusters (alkanes). 
In practice \xn is uniformly sampled in this chemical potential range by step size of 0.05 
eV. The mating operations performed with parents of different number of hydrogen atoms 



may provide a superior sampling of the potential energy landscape 14]. Each GA run has 
typically 10000 generations with 15 random pairs of molecules mated for each generation. At 
the end of each GA run, the candidate structures are further relaxed by ah initio quantum 
chemistry method. 

Recently it has been pointed out that density functional theory (DFT) is unreliable 
for computing hydrocarbon isomer energy differences jlol 16]. Fig JTJ shows the correlation 
between experimental isomer energy differences and those calculated at second order M0ller- 
Plesset perturbation theory (MP2), DFT-B3LYP, DFT-PBEOP levels with basis set of 6- 
31(d) for 9 C 6 H 10 , 7 C 6 H 12 , 5 C 6 H U , 4 C 7 H U and 9 C 7 H 16 molecules. MP2 gives a best 
correlation with the experiments, especially for the low-energy molecules. Hence all the 
candidate hydrocarbon structures in the final pool from GA are relaxed at level of MP2 / 6- 



31G(d) using GAMESS package[17|. 



RESULTS 



We systematically searched for the lowest energy structures of hydrocarbon molecules 
C m H n (m=l-18; n=0-2m+2) in a wide hydrogen chemical potential range of (-6eV, 0). 
Figj2llT4lshow MP2- relaxed lowest-energy structures from our unbiased global searches. Sev- 
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FIG. 1: (Color online) Computed isomer energy differences at MP2 (black square), DFT-B3LYP 
(red circle) and DFT-PBEOP (blue triangle) level versus experimental values for 34 hydrocarbon 
molecules. 



eral low-energy isomers were also shown for comparison. The ground state geometries of 



pure carbon clusters have been extensively studied and reviewed |1 81] . It usually requires 
higher level of quantum chemistry calculations to determine the fine structures and energy 
differences between the isomers. In this paper we focus on providing the lowest-energy 
hydrocarbon clusters in various hydrogen chemical potentials and the lowest-energy pure 
carbon clusters are intended to be chosen as reference systems. Nevertheless, the general 



features (i.e., chain or ring structure) of the grounc 



are consistent with the results in the literature 



state geometries of pure carbon clusters 



19] . 



Fig. [2] shows the lowest energy structures with a backbone of 6 carbon atoms in the full 
hydrogen chemical potential range. The relative energy E associated with each structure 
is defined as 

E = E CnHm - E Cn - mE H (2) 

where Ec n n m , Ec n and Eh are the total energies of the hydrocarbon C n H m , pure car- 
bon cluster C n and hydrogen atom H, with Eh = — 1.356eV. The energies are eval- 
uated at MP2/6-31G(d) level. The pure C 6 cluster is chosen as the reference point. 
In this C6 group, the lowest-energy hydrocarbon structures in the hydrogen chemical 
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FIG. 2: 11 lowest energy structures of C§H m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV.ll lowest energy structures of 
C§H m in full hydrogen chemical potential range. The number below each structure is the relative 
energy Eo in unit of eV. 

potential range of (-G.OOOeV, -4.155eV), (-4.155eV, -4.073ey), (-4.073eV, -2A73eV), 
(-2.473eV,-2.363eV), (-2.363eV,0) are C 6 (chain), C 6 H 4 (Benzyne), C 6 H 6 (Benzene), 
CqHi2 (a,Cyclohexane) and CqH^ (a,2,2-Dimethylbutane), respectively. Several low energy 
isomers in C§H\2 and CqHu group are also shown for reference. Cyclohexane with boat 
conformation CqH 12(c) is 0.287eV^ higher in energy than Cyclohexane with chair conforma- 



tion C$H 12(3), close to the experimental result 0.238 eV[20j. Methylcyclopentane C^H^ih) 
is in between and 0.197 eV higher than CqH\2{&), close to experimental value 0.177 eV 21], 
CqH 12(d) is rather high in energy because of its diradical structure. The energy orders of the 
four alkane isomers are correctly predicted by MP2 as compared with experimental results 



although the energy difference can be quite small 



2l| 
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FIG. 3: 7 lowest energy structures of C 7 H m in full hydrogen chemical potential range. The number 
below each structure is the relative energy Eo in unit of eV. 

Fig. [3] shows the lowest energy structures with a backbone of 7 carbon atoms in the full 
hydrogen chemical potential range. In this C7 group, the lowest-energy hydrocarbon struc- 
tures in hydrogen chemical potential range of (— 6.000eV, — 3.758eV), (— 3.758eV, — 2A67eV), 
(-2.467eV,-2.314eV), (-2.314eV, 0) are C 7 (chain), C 7 H 8 (Methylbenzene), C 7 H U 
(Methylcyclohexane) and C 7 Hiq (a,2,2-Dimethylpentane), respectively. From C6 to C7 
group, the system mainly choose to add one methyl to the dominant lowest energy hydro- 
carbon molecules (i.e., Benzene and Cyclohexane) and exhibit four different lowest energy 
chemical compositions in the whole hydrogen chemical potential range. The alkane subgroup 
is different and still the one with one end methyl pair has lowest energy, in consistence with 



experiment 



2l[. 



CsH m 

Fig. H] shows the lowest energy structures with a backbone of 8 carbon atoms 
in the full hydrogen chemical potential range. In this C8 group, the lowest-energy 
hydrocarbon structures in hydrogen chemical potential range (— 6.000eV, — 4.243eV), 
(-4.243eV,-2.905eV), (-2.905eV, -2.459eV), (-2.459eV, -2.266), (-2.266eV,0) are C 8 
(ring), C$H 6 (a,Ethynylbenzene), C%H W (a,l,2-Dimethylbenzene), C$Hiq (a,cis-l,3- 
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FIG. 4: 15 lowest energy structures of C%H m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

Dimethylcyclohexane) and CgHi S (a,2,2-Dimethyhexane), respectively. Ethynylbenzene has 
lower energy than Cgifg^Pentalene) because of its conjugated bond configuration. The 
first three low energy Benzene-based isomers in C8H10 subgroup are very close in energy. 
MP2-level calculations can not correctly predict the energy orders between them and fa- 
vor 1,2-Dimethylbenzene. In contrast, gas phase thermochemistry measurement shows that 
CgHio (c,l,3-Dimethylbenzene) is 0.018 eV lower in energy than 1,2-Dimethylbenzene 21]. 



On the other hand, MP2 predicts that C$H W (d,Ethylbenzene) is higher in energy than the 



Dimethylbenzenes in agreement with experiment [2 1|. For the C8H16 subgroup, MP2 also 
predicts that Dimethylcyclohexanes have lower energies than Ethylcyclohexane and the en- 
ergy differences between Dimethylcyclohexanes can be very small. In the alkane subgroup, 
the one with one end methyl pair is still the lowest energy isomer. 
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FIG. 5: 15 lowest energy structures of CgH m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

Fig. [5] shows the lowest energy structures with a backbone of 9 carbon atoms in the full 
hydrogen chemical potential range. In this C9 group, the lowest-energy hydrocarbon struc- 
tures in hydrogen chemical potential range (— 6.000eV, — 3.981eV), (— 3.981eV, — 3.397eV), 
(-3.397eV,-2.563eV), (-2.563eV, -2.454eV), (-2.454eV, -2.240eV), (-2.20eV,0) are 
C 9 (chain), C 9 H 7 , C 9 H 8 (indene), C 9 Hi 2 (a,l,2,5-Trimethylbenzene), C 9 H 1S (1,3,5- 
Trimethylcyclohexane) and Cgi^o (a,2,2,5-Trimethyhexane), respectively. The group shows 
more ground state chemical compositions in the full hydrogen chemical potential range. 
MP2 predicts that the isomer with 1,3,5-Trimethyl has the lowest energy in both C9H12 
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FIG. 6: 7 lowest energy structures of CioH m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

and C9H18 subgroups due to maximal reduction of the repulsion energy between the hydro- 
gen atoms. The ground state geometry of the alkane subgroup is predicted to have three 
methyls in accordance with experiment [211] . 



CioH rn 

Fig. [6] shows the lowest energy structures with a backbone of 10 carbon atoms in the full 
hydrogen chemical potential range. In this C10 group, the lowest-energy hydrocarbon struc- 
tures in hydrogen chemical potential range (— 6.000eV, — 3.868eV), (— 3.868eV, — 2A9SeV), 
(-2.498eV,-2.303eV), (-2.303ey, -2.285eK), (-2.285eV, 0) are C 10 (ring), C W H 8 (Naph- 
thalene), CiqHxq (Adamantane), CiqH 2 o (a,trans-l,2,4,5-Tetramethylcyclohexane) and 
C10H22 (a,2,2,5,5-Tetramethyhexane), respectively. Interestingly, the smallest 3D fragment 
of diamond (i.e., Adamantane) appears as a ground state geometry in the group. The cor- 
responding chemical potential range is quite narrow though. The 2D cyclohexane-based 
structure is still present. The isomer with methy pair at both ends is predicted to have the 
ground state geometry in the alkane subgroup. 
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FIG. 7: 11 lowest energy structures of C\\H m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

CnH m 

Fig. [7] shows the lowest energy structures with a backbone of 11 carbon atoms 
in the full hydrogen chemical potential range. In this Cll group, the lowest-energy 
hydrocarbon structures in hydrogen chemical potential range (— 6.000eV, — 4.0496^), 
(-4.049eV,-3.288eV), (-3.288eV, -2.504eF), (-2.504eV, -2.264eV), (-2.264eV,0) are 
Cu ( rm g), CnH 7 , CuH w (a,2-Methylnaphthalene), CnHis (a,l-Methyladamantane) and 
C11H24 (a,2,2,5,5-Tetramethyheptane), respectively. The energy order of the isomers in 
C11H10 subgroup is correctly predicted by MP2 compared with experiment although the 
energy difference is tiny |2l| . MP2 favors the growth of one methyl at the monohydrogen 
site and predicts that 1-Methyladamantane is 0.169 eV lower in energy than CuHi 8 (b,2- 
Methyladamantane), close to the experimental value of 0.206 eV 21]. One may also observe 
that the 2D cyclohexane-based structure disappears in the ground state geometries of the 
hydrocarbons with carbon number larger than ten. Thus we observe a crossover between 
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FIG. 8: 13 lowest energy structures of C\2H m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

populations of 2D cyclohexane-based structure and 3D diamond fragment with increasing 
size of hydrocarbon. The ground state isomer in the alkane subgroup is still predicted to 
the one with methyl pair at both ends. 



Fig. |8] shows the lowest energy structures with a backbone of 12 carbon atoms in the full 
hydrogen chemical potential range. In this C12 group, the lowest-energy hydrocarbon struc- 
tures in hydrogen chemical potential range (— 6.000eV, — 4.198eV), (— 4.198eV, — 2.561eV), 
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FIG. 9: 10 lowest energy structures of C\^H m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV.10 lowest energy structures of 
C\zH m in full hydrogen chemical potential range. The number below each structure is the relative 
energy Eo in unit of eV. 

(-2.561eV,-2.511eV), (-2.511eV, -2.230eV), (-2.230eV,0) are C 12 (ring), C 12 H 8 (Ace- 
naphthylene) , C12H12 (a,2,6-Dimethylnaphthalene), C12H20 (a,l,3-Dimethyladamantane) 
and Ci2-f^26 (a,2,2,7,7-Tetramethyloctane), respectively. Similar rules as mentioned above 
can be applied here: 2,6-Dimethylnaphthalene is the lowest energy isomer because of the 
maximal reduction of the repulsion energy between the hydrogen atoms; monohydrogen site 
is the preferred site to grow methyl on adamantane; the isomer with methyl pair at both 
ends is still lowest in energy in the alkane subgroup. 
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Cy$H m 



Fig. |9] shows the lowest energy structures with a backbone of 13 carbon atoms in the full 
hydrogen chemical potential range. In this C13 group, the lowest-energy hydrocarbon struc- 
tures in hydrogen chemical potential range (— 6.000eV, — 4.025eK), (— 4.025eV, — 2.516eV), 
(-2.516eV,-2.195eV), (-2.195eV, 0) are C 13 (ring), C 13 H 9 (Phenalene), C 13 H 22 (a,l,3,5- 
Trimethylnaphthalene) and Ci 3 H 2 s (a,2,2,8,8-Tetramethylnonane), respectively. Phenalene 
is the lowest energy isomer in C13H9 subgroup due to its closely packed aromatic rings. 1,3,5- 
Trimethylnaphthalene is predicted to be much lower in energy than C\ 3 H 2 2 (b,dodecahydro- 
lH-Phenalene). The ground state structure in the alkane subgroup still keeps the single 
carbon chain with methyl pair at both ends. 

Fig. [10] shows the lowest energy structures with a backbone of 14 carbon atoms in the full 
hydrogen chemical potential range. In this C14 group, the lowest-energy hydrocarbon struc- 
tures in hydrogen chemical potential range of (— 6.000eV, — 4.113eV), (— 4.113eV, — 3.283eV^), 
(-3.283eV, -2.515e1/), (-2.515eV, -2.399eV), (-2.399eV, -2.16SeV), (-2.163eV,0) are 
Ci4 (ring), CliHg (Cyclopent[fg]acenaphthylene), C U H 10 (a,Phenanthrene), C u H 2 o 
(a,Diamantane), C\±H24 (a,l,3,5,7-Tetramethyladamantane) and C U H 30 (a, 2, 2, 5,8,8- 
penmethylnonane), respectively. Contrary to the rule of maximal reduction of repulsion 
energy between hydrogen atoms, Phenanthrene is correctly predicted to be 0.283 eV lower 
in energy than Ci4-£/io(b,Anthracene), comparable to the experimental result of 0.216 eV. 
A bigger diamond fragment (Diamantane) appears as ground state structure in a narrow 
hydrogen chemical potential range. The fifth methyl appears in the middle of the carbon 
chain besides methyl pair at both ends of the lowest energy isomer in the alkane subgroup. 

Fig. [UJ shows the lowest energy structures with a backbone of 15 carbon atoms 
in the full hydrogen chemical potential range. In this C15 group, the lowest-energy 
hydrocarbon structures in hydrogen chemical potential range of (— 6.000eV, — 4.239eV^), 
(-4.239eV,-2.729ey), (-2.720eV, -2.519eV), (-2.519eV, -2.243eV), (-2.243eV,0) are 
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FIG. 10: 12 lowest energy structures of Cy^H m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

C15 (ring), C15H9 (a), C\%Hvi (3-Methylphenanthrene), C15H22 (a,3-Methyldiadamantane) 
and Ci 5 H 32 (a,2,2,5,5,8,8-Hexamethylnonane), respectively. Ci 5 H 9 (a) is predicted to be the 
lowest energy isomer due to its compact structure and minimal influence of the five fold ring. 
Methyl is still preferred to grow at the monohydrogen site of diadamantane. Three pairs of 
methyls are distributed symmetrically on the carbon chain in the ground state structure of 
alkane subgroup. 
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0.000 
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-37.135 -37.001 
(c) (d) 
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FIG. 11: 10 lowest energy structures of C\^H m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

C16 C16H10 






0.000 (a) -40.507 (b) -38.975 
C16H24 C16H34 




-74.928 (a)-97.153 (b)-97.079 



FIG. 12: 6 lowest energy structures of C\§tL m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

Fig. [T2l shows the lowest energy structures with a backbone of 16 carbon atoms in the full 
hydrogen chemical potential range. In this C16 group, the lowest-energy hydrocarbon struc- 
tures in hydrogen chemical potential range (— 6.000eV, — 4.051eV), (— 4.051eV, — 2.459eV), 
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C17H26 C17H36 
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-82.376 (a)-104.501 (b)-104.490 



FIG. 13: 6 lowest energy structures of CnH m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eq in unit of eV. 



(-2.459eV, -2.223eV), (-2.223eV, 0) are C 16 (ring), C W H W (a,Pyrene), C 16 H 2 4, and C 16 H M 
(a, 2, 2, 5, 5, 9, 9-Hexamet hylo ct ane ) , resp ect ively. 



Fig. [TBI shows the lowest energy structures with a backbone of 17 carbon atoms in the full 
hydrogen chemical potential range. In this C17 group, the lowest-energy hydrocarbon struc- 
tures in hydrogen chemical potential range of (— 6.000eV, — 4.046ey), (— 4.046eV, — 3.496eV^), 
(-3.496eV, -2.456eV), (-2.456eV, -2.213eV), (-2.213eV, 0) are C 17 (ring), C 17 H U , C 17 H 12 
(2-Methylpyrene), Ci7#26> and CyjH^Q (a,2,2,6,6,10,10-Hexamethyludecane), respectively. 



Fig. shows the lowest energy structures with a backbone of 18 carbon atoms 
in the full hydrogen chemical potential range. In this C18 group, the lowest-energy 
hydrocarbon structures in hydrogen chemical potential range (— 6. OOOeV, —4.377eV), 
(-4.377ey,-2.564eK), (-2.564eV, -2.506eV), (-2.506eV, -2.234eF), (-2.234eV,0) are 
Cig (ring), Ci$H w (a), Ci 8 H u (2,7-Dimethylpyrene), Ci 8 if 2 4 (Diamantane) and Ci$H 38 
(a,2,2,5,5,ll,ll-Hexamethyldodecane), respectively. Similar empirical rules are observed 
in C16-C18 groups: Poly-aromatic rings tend to form compact structure; Methyl favors 



C17H, 



rn 



in 
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0.000 (a)-43.769(b)-43.606 -54.02( 
C18H24 C18H38 






-79.087 (a)-l 10.361 (b)-l 10.284 



FIG. 14: 7 lowest energy structures of C\%H m in full hydrogen chemical potential range. The 
number below each structure is the relative energy Eo in unit of eV. 

the monohydrogen site in diamond fragment structure; ground state structure of alkane 
subgroup favors methyl pairs spaced by at least two carbon atoms over the main carbon 
chain. 



CONCLUSION 



An unbiased evolution-based optimization method combined with Brenner's empirical 
potential is used to search for ground state structures of hydrocarbon molecules in a wide 
hydrogen chemical potential range. The resultant structures are further sorted by quantum 
chemical calculations at MP2 level. The collection of lowest energy structures of hydrocarbon 
molecules C m H n (m=l-18; n=0-2m+2) is presented. A crossover between populations of 2D 
cyclohexane-based structure and 3D diamond fragment with increasing number of carbon 
atoms of the hydrocarbon molecules is demonstrated. Besides the PAH compounds, we 
also show that PAH with methyls can also be important in the interstellar medium. The 
spectra of the PAH compounds and how are they affected by the various configurations of 
the radicals (e.g., methyl) are worthy of further studying. 
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